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Light-cone path integral approach to the induced radiation in QED and 
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I discuss the basic ideas of the light-cone path integral approach to the induced radiation in QED and QCD 
and recent applications to the induced parton energy loss. 



1. The induced radiative energy loss and the 
Landau-Pomeranchuk-Migdal (LPM) effect [Tl2] 
in QED and QCD have attracted much attention 
in recent years, see |3I4I5| and references therein. 
It is mainly because of the first accurate measure- 
ments of the LPM effect in QED at SLAC 0, and 
the possibility to use jet quenching for probbing 
the quark-gluon plasma (QGP) produced in high- 
energy ^.^-collisions |4I7I8I9| . The most general 
approach to the induced radiation applicable in 
both QED and QCD is the light-cone path inte- 
gral (LCPI) approach [TUJ (see also |11I12I4| L It 
accurately treats the mass and finite-size effects, 
and applies at arbitrary strength of the LPM ef- 
fect. This approach gives excellent description of 
the SLAC and SPS [J3] data on the photon 
radiation from high-energy electrons [JJ| . In this 
talk I discuss the basic concepts and recent appli- 
cations of the LCPI approach. 

2. The starting point of the LCPI formalism is 
the representation of the wave functions of high- 
energy free particles in the form (it is assumed 
that the angles between the z axis and velocities 
are small) 

V>i(r) = cxpiiEjZpj^ip, z) , (1) 

where p — r_L, Uj is a spin operator, and <pj sat- 
isfies the Schrodinger equation 
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Figure 1. 



The Schrodinger mass pj in J3J) equals the parti- 
cle energy Ej . The z-evalution of 4>j is described 



by the Green's function for the Schrodinger equa- 
tion lf2"|. Kj. In vacuum the amplitude squared 
|(6c|T|a)| 2 for the a — > be transition is given 
by the diagram shown in Fig. la, where the — > 
(<— ) lines correspond to Kj (Kj), the dashed 
lines show the initial (at Zi ~ — oo) and final (at 
Zf ~ oo) transverse density matrices, pi, pf, and 
the integration over the transverse coordinates of 
the endpoints (and vertices) and the longitudinal 
coordinates of the vertices (below z\ and zi for 
upper and lower parts) is implicit. In vacuum 
(bc\T\a) = (if the particle a is not produced in 
a hard reaction, and ui a ^ mi) -f-77i c ), however, in 
an external potential it is not the case. 

To evaluate the induced a — > be transition in 
an amorphous medium one should sum i-channel 
exchanges between the fast particles and parti- 
cles in the medium and perform averaging over 
the medium states. The key idea of the LCPI 
approach is to represent all Kj in the Feynman 
path integral form, and perform averaging over 
the potential at the level of the integrands be- 
fore integrating over the trajectories. After av- 
eraging over the medium states the interaction 
of the fast particles with the medium turns out 
to be translated to an effective interaction be- 
tween trajectories. Both in QED and QCD (note 
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that in QCD we work at the level of two-gluon 
i-channel exchanges) from the viewpoint of this 
interaction the trajectories corresponding to K* 
may be viewd as antiparticle trajectories (we call 
— > and <— lines "particle" and "antiparticle" tra- 
jectories). In medium the effective Lagrangian for 
the path integral form of the diagram la reads 



P'a,P'a,P', z '\S\Pa,Pa,Pi z ) = K a (p' a , APai Z ) 



Leff — L q (t p ) — L {fp) + L int (T p , Tp), 



(3) 



where t p , Tp are the sets of the transverse co- 
ordinates for "particles" and "antiparticles" , Lq 
and Lq are the corresponding free Lagrangians, 
say, Lq(t p ) — fi-iPi /2- The interaction term 
reads L int (T p ,T p ) = in(z)a x (Tp,T p ) /2, where 
ax is the diffractive operator for X = "parti- 
cles" +" antiparticles" system scattering off a par- 
ticle in medium, n(z) is the number density of the 
medium. In QED ax is simply the cross section, 
but in QCD for the 4-body part [z 2 < z < Zf) 
ax acts as an operator in color space. The one- 
particle spectrum, say, integrated over q c (here- 
after q stands for the transverse momentum), 
may be represented in a form which does not con- 
tain the 4-body part at all. Since integration over 
q c gives 5{p c — p g ), the diagram la is transformed 
into lb in this case. The spectrum correponding 
to the diagram lb reads 
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X g(Pf\S b - b ® S bca ® S a a\p t ), (4) 

where x = x b — Pb,z/Pa,z, Q is the vertex factor, 
Saa, S bb , S bc5 are the evolution operators for the 
corresponding L e ff. For the two-body parts the 
path integrals can be taken analytically, say, 

(Pa,p'a,z'\Saa\P a ,Pa> z ) = K a{p' a , z '\p a i z) 
xK *a(p'a' Z '\Pa, z )®aa(p' a ~ p'a^'\Pa ~Pai Z )i ( 5 ) 

where $ aa =exp [-■ | f*dzn(z)a a a(p a ( z ) - Pa( z ))} 
should be evaluated for the straight trajectories 
of a and a. For the 6ca-part one can write 
/ Dp b Dp c Dp a as / DpDp a Dp a , where p = p b 
p c , and p a = x b p b + x c p c is the center-of-mass 
coordinate of the be system. The p a and p a in- 
tegrals can be taken analytically, and in the new 
variables S bc a can be represented as 



xK a(p'a, z '\Pai *)K{P\ Z '\P, z ), 



(6) 



where K, is the Green's function for the Hamilto- 
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H =-m*)[dp-) +«*>r)+r t - (7) 

Here v(z,p) = -in(z)a bca {p, p a - p a )/2 should 
be evaluated for the straight trajectories p a , p a , 
M(x) = E a x(l-x), L f = 2M(x)/e 2 is the forma- 
tion length, e 2 = [m b x c + m^x b — m^XbX c ]. Hav- 
ing JSJ and ©, one can obtain 

^ b = ^^jd^M-i^ ■ ^J^ W 

z f 

\ I dz 2 g<Pf(T,z 2 )lC(T,z2\0,z 1 )<S> i (xT,z 1 ), (8) 
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The contribution stemming from the region of 
large \zi^\ in © can be expressed via the light- 
cone wave-function of the a — ► be transition in 
vacuum, \& a c , 0]. Then JSJl takes the form 

dP 



^^Re I dr exp(-iq 6 r) 



dxdq b (27r) 2 

r z f r z f 

x dzi / dz 2 g{® f {T,z 2 )[}C(T,z 2 \0, zi) 

J Zi •* Z\ 

-K, v (T,Z 2 \0,Zl)]^i(xT,Z 1 ) 
+ [<f>f(T,Z 2 ) - l]K v (T,Z 2 \0,Z 1 )[^i{xT,Z 1 ) - 1]} 

+7^)2 / drdr' exp(-iq b r)^ b a c *(x,T' - t) 
x * c (ar,r / )[$/(r,z i ) + $ i (xT,^)-2] 



(9) 



where IC V is the Green's function for v = 0. For 
Lf ^> L (L is the thickness of the medium) using 
the representation of ^ a c via the z-integral of JC 
esteblished in one can represent J3J as 



-^-rRe/ drdr' exp(— iq b r) 

iny J 



dP 

dxdq b (2t 

^*{x,t'-t)^{x,t') 

x [2r bca (x', XT)-r bb {r)-r a - a {xr)} , (10) 
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where T h = 1 - exp [-^ dzn(z)}. 

From © one obtains for the x-spectrum ^U] 

dP f Zf f Zf 

— = 2ReJ dzij dz 2 g [K(p 2 , z 2 \p 1} Zi) 

- K v (p 2 ,z 2 \p 1 ,zx)] ■ (11) 

Pi=P 2 =0 

In [15] . separating the N = 1 rescattering, we 
have represented Ijlljl as a sum of the Bethe- 
Hcitler spectrum plus an absorptive correction re- 
sponsible for the LPM effect. This form has been 
used ^3] for successful description of the SLAC 
[Hj and SPS Q3] data on the LPM effect in photon 
bremsstrahlung from high energy electrons. 

For a particle produced in the medium it is con- 
venient to rewrite Qllf) in another form. For gluon 
emission from a quark this new form reads |16j 
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Figure 2. The quark energy loss in the nuclear 
matter (a) for L — 3, 6, 9 fm and expanding QGP 
(b) for L = 6 fm for RHIC and LHC. 



da»/}(x,z) 



Re / dz\ I dz 2 I dpg 



dx j 
xlC v (p 2 , z 2 \p, z)a 3 (p)K,(p, z\p 1 ,z 1 ) 



. (13) 



where 03 = <7 gq q, L is the quark pathlength in 
the medium. The daffi/dx (O can be viewed 
as an effective Bethe-Heitler cross section which 
accounts for the LPM and finite-size effects. Ne- 
glecting spin-flip transition Ijl3(l can be written as 
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where P(x) 9q 



Cf[1 + (1 — x) 2 ]/x is the usual 
splitting function, and F is the solution to the ra- 
dial Schrodinger equation for the azimuthal quan- 
tum number m = 1 
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The boundary condition for F(£,p) reads F(£ = 
0, p) = y/p<73(p)eKi(ep), where K\ is the Bessel 
function. 

3. The formulas l|12f> . I|14() are convenient 
for numerical calculations with an accurate 
parametrization of 03. Note that the widely 
used oscillator approximation 03 (p) cx p is too 
crude for the induced gluon emission [17116) . We 
perform calculations for running a s frozen at 
a s — 0.7 at low momenta (for incorporation of 
running a s see |lfijL In Fig. 2 we show the en- 
ergy dependence of the quark energy loss AE = 
E J dxxdP/dx for the cold nuclear matter and 
for expanding QGP for RHIC and LHC condi- 
tions. For nuclear matter we take m g = 0.75 
GeV obtained from the analysis of the low- a; pro- 
ton structure function within the dipole BFKL 
equation |18ll9j . It agrees well with the inverse 
gluon correlation radius in the QCD vacuum |20] . 
For QGP we use m g — 0.4 GeV obtained in the 
quasiparticle model from the lattice data in |21| . 
For the light quark we take m q = 0.3 GeV. For 
QGP this quark mass was obtained in Note 
that our results are not very sensitive to the light 
quark mass. For the QGP we also show the re- 
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Figure 3. The nuclear modification factor for 
central Au + Au collisions at y/s = 200 GeV for 
quark (solid line) and gluon (dashed line) jets ob- 
tained with m g = 0.4 GeV (thick lines) and m g = 
0.75 GeV (thin lines). The experimental points 
(from (2H|) are for: circle - Au + Au — > tt° + X 
(0-10% central) [PHENIX Collaboration], square 
- Au + Au -> h ± + X (0-10% central) [PHENIX 
Collaboration], star - Au + Au -> h ± + X (0-5% 
central) [STAR Collaboration]. 



suits for c-quark, m q = 1.5 GeV. The Debye 
screening mass, /id, in QGP was obtained from 
the perturbative relation po = V2m g « 0.57 
GeV. We use the Bjorken model for the QGP 
expansion with T 3 r = TqTq and the initial con- 
ditions suggested in |23]: T = 446 MeV and 
r = 0.147 fm for RHIC (for Au + Au at y/1 = 200 
GeV), and T = 897 MeV and r = 0.073 fm for 
LHC (Pb + Pb at ^ = 5.5 TeV). For the QGP 
we take L = 6 fm. It is a about the life-time of 
the QGP (and mixed) phase r max ~ Ra for cen- 
tral heavy-ion collisions, since due to the trans- 
verse expansion the hot QCD matter should cool 
quickly at t > Ra ■ 

AE for nuclear matter obtained using H12fl - I|15|) 
with an accurate parametrization of the imagi- 
nary potential by a factor of about 2 smaller than 
that obtained previously [H] in the oscillator ap- 
proximation. Our AE for nuclear matter is con- 
siderably smaller than prediction of [21] in the 
HERMES [23 energy region E < 20 GeV. The 
small AE shows that for the HERMES condi- 
tions the effects of hadron absorption and string 
tension should be more important than the gluon 
emission since the hadron formation time If is 
about the nuclear size. Note that the estimate 



If ~ 40 fm for E = 8 GeV given in [2U is abso- 
lutely unrealistic. The authors of [21] do not pay 
any attention to the fact that a quark with such 
an energy should be stopped at L < 8 fm since 
the string tension is ~ 1 GeV/fm. 

One can see from Fig. 2b that for LHC condi- 
tions the jets with E < 20 GeV should practically 
be absorbed in the QGP. It means that the sur- 
face jet production dominates. The energy loss 
for c-quark is smaller than that for light quarks 
by a factor of ~ 2 at E ~ 5 - 10 GeV. 

The effect of parton energy loss on the high-pT 
hadron production in A+A collisions can approx- 
imately be described in terms of effective hard 
partonic cross sections which account for the in- 
duced gluon emission [SJ. Using the power- low 
parametrization for cross section of quark pro- 
duction in p + p collisions cx l/p!^ PT ^ the nuclear 
modification factor 
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can be written as IT? 



Raa{pt)~Pq{ptY 



J(pt)Jq 



dz 



y n{PT)-2j,h, PT 



1 dl(x, -tF-^) 

dx 



(17) 



where J(p T ) = J* dzz n ^- 2 D h q {z, ^f), P is the 
probability of quark propagation without induced 
gluon emission, dI{x,pT)/dx is the probability 
distribution in the quark energy loss for a quark 
with E — pt : Dq(z,PT/z) is the quark fragmen- 
tation function. A convenient parametrization for 
Po and dl(x,px) / dx in terms of dP/dx (similar 
to that used for photon emission in J3]) is given 
in PS]- In Fig. 3 we compare the theoretical Raa 
with the RHIC data on Au + Au collisions at 
■s/s = 200 GeV [2E1- The theoretical curves have 
been obtained for L — 4.9 fm. It is the typical 
parton pathlength in the QGP (and mixed) phase 
for T m ax = 6 fm. The results for the quark and 
gluon jets are shown separately (for y/s = 200 
GeV the quark and gluon contributions are com- 
parable). The suppression is somewhat stronger 
for gluon jets. As for Fig. 2b the calculations are 
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Figure 4. The kinematic 7Vfactor for N = 1 
rescattering for the QGP, L = 5 fm, n(z) = const. 
Solid line shows our results. The points show the 
GLV predictions. 



performed for m q — 0.3 GeV and m g = 0.4 GeV. 
To illustrate the m g -dependence we also show the 
results for m g = 0.75 GeV (with the same \xd as 
for m g = 0.4 GeV), which seems to be reason- 
able for the mixed phase and for gluons with large 
formation length Lf > L. The above two values 
of m g give reasonable lower and upper limits of 
the infrared cutoff for the induced gluon emission 
for RHIC conditions. The theoretical Raa for 
m g = 0.4 GeV is in quite good agreement with 
the experimental one. 

4. The LCPI approach neglects the kinematic 
bounds. The GLV group reported that the 
finite kinematic limits suppress strongly the en- 
ergy loss even at E ~ 500 GeV. To check it in 
[2*2*] we have calculated the dominating N = 1 
rescattering contribution to the quark energy loss 
using the ordinary diagram language. We have 
obtained a small kinematic effect. To demon- 
strate it in Fig. 4 we show the kinematic if -factor 
K = AEf.i./AEu., where AE f . L and AJSj.j. are 
the energy losses for finite and infinite kinematic 
limits. Our calculations are performed for the 
kinematic limits as in [""", also similarly to [Jj 
we use fixed ct s , n(z) = const, L = 5 fm, and 
Tn g = (J-d = 0.5 GeV. One sees that our if -factor 
contrary to the GLV predictions is close to unity 
even for E ~ 5 GeV. It says that the LCPI ap- 
proach has a quite good accuracy for E > 5 GeV. 

I thank the Organizers of Difftaction'04 for the hos- 
pitality and financial support during the Workshop. 



REFERENCES 

1. L.D. Landau and I.Ya. Pomeranchuk, Dokl. 
Akad. Nauk SSSR 92 (1953) 535, 735. 

2. A.B. Migdal, Phys. Rev. 103 (1956) 1811. 

3. S.R. Klein, Rev. Mod. Phys. 71 (1999) 1501. 

4. R. Baier, D. Schiff and B.G. Zakharov, Ann. 
Rev. Nucl. Part. 50 (2000) 37. 

5. A. Kovner and U.A. Wiedemann, 
|hep-ph/0304151| 

6. P.L. Anthony, R. Becker-Szendy, P.E. Bosted, 
et ah, Phys. Rev. D56 (1997) 1373. 

7. M. Gyulassy, P. Levai and I. Vitev, Nucl. 
Phys. B594 (2001) 371. 

8. R. Baier, Yu.L. Dokshitzer, A.H. Mueller and 
D. Schiff, JHEP 0109 (2001) 033. 

9. B.G. Zakharov, JETP Lett. 65 (1997) 615. 

10. B.G. Zakharov, JETP Lett. 63 (1996) 952. 

f f . B.G. Zakharov, Phys. Atom. Nucl. 61 (1998) 
838. 

12. B.G. Zakharov, JETP Lett. 70 (1999) 176. 

13. H.D. Hansen, U.I. Uggerh0j, C. Biino et al., 
Phys. Rev. Lett. 91 (2003) 014801. 

14. B.G. Zakharov, Phys. Atom. Nucl. 62 (1999) 
1008; JETP Lett. 78 (2003) 759. 

15. B.G. Zakharov, JETP Lett. 64 (1996) 781. 

16. B.G. Zakharov JETP Lett. 80 (2004) 617; 
|hep-ph/0410321| 

17. B.G. Zakharov, JETP Lett. 73 (2001) 49. 

18. N.N. Nikolaev, B.G. Zakharov and 

V.R. Zoller, Phys. Lett. B328 (1994) 486. 

19. N.N. Nikolaev and B.G. Zakharov, Phys. 
Lett. B327 (1994) 149. 

20. E.V. Shuryak, Rev. Mod. Phys. 65 (1993) 1. 

21. P. Levai and U. Heinz, Phys. Rev. C57 (1998) 
1879. 

22. B.G. Zakharov, JETP Lett. 80 (2004) 67. 

23. R.J. Fries, B. Miiller and D.K. Srivastava, 
Phys. Rev. Lett. 90 (2003) 132301. 

24. E. Wang and X.N. Wang, Phys. Rev. Lett. 89 
(2002) 162301. 

25. A. Airapetian et al. (HERMES Collabora- 
tion), Phys. Lett. B577 (2003) 37. 

26. D. d'Enterria, Invited overview talk at 39th 
Rencontres de Moriond on QCD and High- 
Energy Hadronic Interactions, La Thuilc, 
Italy, 28 Mar - 4 Apr 2004; nucl-ex/0406012 



